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Further Properties of Wireless Channel Capacity 

Fengyou Sun and Yuming Jiang 


Abstract 

Future wireless communication calls for exploration of more efficient use of wireless channel capacity to meet the increasing 
demand on higher data rate and less latency. However, while the ergodic capacity and instantaneous capacity of a wireless 
channel have been extensively studied, they are in many cases not sufficient for use in assessing if data transmission over the 
channel meets the quality of service (QoS) requirements. To address this limitation, we advocate a set of wireless channel 
capacity concepts, namely “cumulative capacity”, “maximum cumulative capacity”, “minimum cumulative capacity”, and “range 
of cumulative capacity”, and for each, study its properties by taking into consideration the impact of the underlying dependence 
structure of the corresponding stochastic process. Specifically, their cumulative distribution function (CDFs) are investigated 
extensively, where copula is adopted to express the dependence structures. Results considering both generic and specific dependence 
structures are derived. In particular, in addition to i.i.d., a specially investigated dependence structure is comonotonicity, i.e, the 
time series of wireless channel capacity are increasing functions of a common random variable. Appealingly, copula can serve 
as a unifying technique for obtaining results under various dependence assumptions, e.g. i.i.d. and Markov dependence, which 
are widely seen in stochastic network calculus. Moreover, some other characterizations of cumulative capacity are also studied, 
including moment generating function, Mellin transform, and stochastic service curve. With these properties, we believe QoS 
assessment of data transmission over the channel can be further performed, e.g. by applying analytical techniques and results of 
the stochastic network calculus theory. 


I. Introduction 

In future wireless communication, there will be a continuing wireless data explosion and an increasing demand on higher 
data rate and less latency. It has been depicted that the amount of IP data handled by wireless networks will exceed 500 
exabytes by 2020, the aggregate data rate and edge rate will increase respectively by 1000 x and 100 x from 4G to 5G, and the 
round-trip latency needs to be less than 1ms in 5G m. Evidently, it becomes more and more crucial to explore the ultimate 
capacity that a wireless channel can provide and to guarantee pluralistic quality of service (QoS) for seamless user experience. 

Information theory provides a framework for studying the performance limits in communication and the most basic measure 
of performance is channel capacity, i.e., the maximum rate of communication for which arbitrarily small error probability can be 
achieved 0 - Due to the time variant nature of a wireless fading channel, its capacity over time is generally a stochastic process. 
To date, wireless channel capacity has mostly been analyzed for its average rate in the asymptotic regime, i.e., ergodic capacity, 
or at one time instant/short time slot, i.e., instantaneous capacity. For instance, the first and second order statistical properties 
of instantaneous capacity have been extensively investigated, e.g. in 0, 0- However, such properties of wireless channel 
capacity are ordinarily not sufficient for use in assessing if data transmission over the channel meets its QoS requirements. 
This calls for studying other properties of wireless channel capacity, which can be more easily used for QoS analysis. To meet 
this need constitutes the objective of this paper. 

Specifically, we advocate in this paper a set of (new) concepts for wireless channel capacity and study their properties. 
These concepts include “cumulative capacity”, “maximum cumulative capacity”, “minimum cumulative capacity”, and “range 
of cumulative capacity”. They respectively refer to the cumulated capacity over a time period, the maximum and the minimum 
of such capacity within this period, and the gap between the maximum and the minimum. 

Among these (new) concepts, the wireless channel cumulative capacity of a period is essentially the amount of data 
transmission service that the wireless channel provides (if there is data for transmission) 0 or is capable of providing 
(if there is no data for transmission) 0 in this period. For the former, the concept is closely related to the (cumulative) service 
process concept that has been widely used in the stochastic network calculus literature, e.g. in |[5l— fT4l . In particular, in these 
works when charactering the cumulative service process using server models of stochastic network calculus and/or applying 
the cumulative service process concept to QoS analysis, some special assumptions on the dependence structure of the process 
are often considered, such as independence 0-0 and Markov property nia, urn, im. 

In addition, we introduce “maximum cumulative capacity”, “minimum cumulative capacity” and “range of cumulative 
capacity” that are new but we believe are also crucial concepts for analyzing QoS performance of wireless channels. This is 
motivated by the fact that, even with the CDF (i.e. full characteristics) or its bounds of the cumulative capacity known, it may 
still be difficult to perform QoS analysis of the channel. (One can easily observe this difficulty by assuming fluid traffic input 
and trying to find backlog bounds from queueing analysis of the channel. See e.g. 0 ). As a special case of these concepts, 
forward-looking and backward-looking variations of them are also defined, which turn out to be useful in different application 
scenarios. 

For the investigation, unlike most existing work in the stochastic network calculus literature, the present paper mainly focuses 
directly on the cumulative distribution functions (CDFs) of the corresponding processes of these (new) concepts. For their other 
characterizations, e.g. moment generating function 0, Mellin transform GU, and stochastic service curve 0, a number of 
results are also reported for cumulative capacity to exemplify how such properties may be analyzed, but this is not focused. An 
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underlying reason is that a random variable is fully characterized by its CDF. To unify the investigation for each concept, we 
introduce copula as a technique to account for the various dependence structures implied by different possible properties of the 
process. For instance, besides a process with i.i.d. increments, a Markov process has also been proved to have a dependence 
structure tm oil- In addition to such dependence assumptions for which many results in the stochastic network calculus are 
available, we use comonotonicity as a dependence structure when there is a strong time dependence in the channel, i.e., when 
the time series of instantaneous wireless channel capacities can be represented as increasing functions of a common random 
variable. Moreover, generic results under arbitrary dependence structures are also obtained when only the marginal distribution 
functions are assumed to be known. 

To remark, the idea of taking advantage of specific dependence structures in analysis can be found in the stochastic network 
calculus literature, e.g., independent increments ( 6 ]|—| 8 ]l and Markov property ITOl . |I2| . Ifl3l . However, such diverse dependence 
structures are investigated separately, without a unified technique. In addition, the literature investigation mainly focuses on 
the stochastic service curve characterization of the cumulative service process and on applying it to QoS performance analysis. 
Little has directly focused on the probability distribution function characteristics of the cumulative capacity. Moreover, to 
the best of our knowledge, there is no previous work focusing on the probabilistic distribution function characteristics of the 
maximum, the minimum and the range of cumulative capacity of a wireless channel. In |18), the concept of copula is brought 
into stochastic network calculus, which is applied to consider the dependence in the superposition property of arrivals. Different 
from m, our focus is the cumulative capacity processes that are related to the cumulative service process, while not the arrival 
process. In addition, we use copula to feature the dependence structures in the considered cumulative capacity processes. 

The contributions of this worlfl are several-fold. (1) Several concepts for studying wireless channel capacity are introduced, 
namely “cumulative capacity”, “maximum cumulative capacity”, “minimum cumulative capacity” and “range of cumulative 
capacity”. For the latter three concepts, we originally introduce their forward-looking and backward-looking versions and 
highlight the fundamental difference between a forward-looking concept and its backward-looking counterpart. (2) A copula 
technique is introduced to unify the analysis under difference dependence conditions. (3) The probability distribution function 
characteristics of the processes corresponding to the introduced concepts are investigated, and various exact solutions or 
bounds on the probability distribution functions are derived. (4) Other characteristics, such as moment generating function, 
Mellin transform and stochastic service curve, of the cumulative capacity process are investigated to exemplify how such 
characteristics of the processes corresponding to the introduced concepts can be studied. Based on these properties, QoS 
assessment of data transmission over a wireless channel can be further performed by, e.g., exploiting the stochastic network 
calculus theory (see e.g. El). 

The remainder of this paper is structured as follows. The fundamental concepts of wireless channel capacity, including 
instantaneous capacity, cumulative capacity, maximum cumulative capacity, minimum cumulative capacity and the range of 
cumulative capacity, are first introduced in Sec. |TT] Also in Sec. [Ill preliminaries for later analysis, including those on copula, 
non-Granger causality and change of measure, are described. In Sec. [Till copula is elaborated as a unifying technique for analysis 
of cumulative capacity processes under different dependence structures. The probability distribution function characteristics, 
particularly the CDF, for cumulative capacity and maximum/minimum cumulative capacity are analyzed respectively in Sec. 
m and Sec. m Other characterizations of cumulative capacity are elaborated in Sec. [VU Finally, the paper is concluded and 
future work is discussed in Sec. [yin 

II. Fundamental Concepts and Preliminaries 

A. Fundamental Concepts 

1) Instantaneous Capacity: Consider a wireless channel. We assume discrete time t = 1,2,..., and that the instantaneous 
capacity lfl9l or mutual information lf20l C\t) of the channel at time t can be expressed as a function of the instantaneous 
SNR 7 1 at this time 0: 

C(t) = log. 2 (g{j t )). (1) 

For single input single output (SISO) channels, if CSI is only known at the receiver, the instantaneous capacity or the mutual 
information of the channel, assuming flat fading, can be expressed as 

C(f)=log 2 (l+ 7 |Mi)| 2 ), (2) 

where h(t) is a stochastic process describing the fading behavior, \h(t)\ denotes the envelope of h(t), 7 = P/NqW denotes 
the average received SNR per complex degree of freedom, P is the average transmission power per complex symbol, Nq/2 is 
the power spectral density of AWGN, and W is the channel bandwidth. 

In the literature, PDF or CDF of the instantaneous capacity is available for various types of channels, e.g. Rayleigh channel 
Ifm Rice channel (22l . Nakagami-m channel 1231 . Suzuki channel (24], and more @. Specifically, the CDF of the Rayleigh 
channel instantaneous capacity is expressed as lf 2 TTl 

F C ( t ){r) = (3) 

lr This is still an on-going work, subject to significant revisions. The contributions are summarized on this on-line version and will be significantly extended 
when more results are added. 
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2) Cumulative Capacity: We define the cumulative capacity through period (s,i] as 

t 

S(s,t) = J2 w, 

2=S + 1 

where C(i) is the instantaneous capacity at time i. 

3) Maximum Cumulative Capacity: We define the maximum cumulative capacity in period (0, i] as 

S{0,t) = sup S(j, k) = sup C i i ) 

where C{i) is the instantaneous capacity at time i. 

Fixing j = 1 in 0, we obtain a forward-looking of the maximum cumulative capacity, i.e., 

3(0,t)= sup 5(0,fc) = sup ( 5^ C'(z) N ] , 

1 <k<t 1 <k<t \~l J 

while fixing k = t in <0, we obtain a backward-looking of the maximum cumulative capacity, i.e., 

*S{Q,t)= sup S(j.t) = sup VC(,] . 

i<f<* i <3<t J 

4) Minimum Cumulative Capacity: We define the minimum cumulative capacity in period (0,t] as 

aM) s = i<j?L (£ c<i) 

\ '-j 

where C(i) is the instantaneous capacity at time i. 

Similarly, fixing j = 1 in ©, we obtain & forward-looking of the minimum cumulative capacity, i.e.. 


S(0,t)= inf 5(0, k) = inf V C(i) 

-> l<k<t \ 


l<k<t 


and fixing k — t in 0 , we obtain a backward-looking of the minimum cumulative capacity, i.e.. 




i <j<t 


. i=3 


5) Range of Cumulative Capacity: We define the range of cumulative capacity in period (0, f] as 

R(0,t) = S(0,t)-S(0,t). 


(4) 


(5) 


( 6 ) 


(7) 


( 8 ) 


(9) 


( 10 ) 


( 11 ) 


The range can also have variations based on the selection of forward-looking and backward-looking expressions of the 
maximum and minimum cumulative capacity, e.g, forward-looking = 1^(0 ,t) — (0,t) and backward-looking 

4(o,i) = ^(o,t)-^(o,t) 

Remark 1. 5(0, t) is essentially a partial sum process in probability theory. When studying the maxima of partial sums 
in probability theory, the forward-looking versions of 5(0, f) and R(0,t) are typically focused 1231. 4261/ . This is probably 
due to that the studied probability problems are often forward-looking in nature, e.g., at time 0 , making decision based on 
possible happenings in the future time t. However, for QoS analysis of communication networks, backw’ard-Iooking is more 
important, e.g. when deciding delay at time t, we have to look backwards to see how much serx’ice the follow has experienced. 
As directly seen from their definitions, the forward-looking maximum cumulative capacity is fundamentally different from its 
backward-looking counterpart. As a consequence, results in probability theory for the process of maxima of partial sums should 
be used with care when they are applied to queueing analysis. In general, under time reversibility assumption, the results for 
the forward-looking definitions may be extended for application to the backward-looking definitions. 

The range describes the gap between the maximum and the minimum, and hence the smaller the range, the closer are the 
two. In particular, when the range is small or approaching 0, the maximum becomes the minimum, hence should also equal 
the cumulative capacity. Range may be used as a measure to characterize the tightness between an upper bound and a lower 
bound on the cumulative capacity. The idea of using a gap between analytical bounds themselves or between the bounds and 
the exact results to characterize the tightness or accuracy of the obtained bounds has been used in e.g. Information Theory to 
investigate / characterize the information capacity of a channel W.I In the context of stochastic network calculus, such an 
idea was exploited in 4251 / to study the accuracy of the obtained bounds. 
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B. Preliminaries 

1) Copula: Copula is a well-known concept for dependence modeling by decoupling the joint distribution function into the 
dependence structure and marginal behavior. 

Definition 1. f29f, g3Ujl A d-dimensional copula is a distribution function on [ 0 ,l] d with standard uniform marginal 

distributions. 


It is equivalent to say that a copula is any function C : [0, l] ra —> [0,1], which has the following three properties: |29l , i30j 

(1) C(ui,..., uf) is increasing in each component Ui. 

( 2 ) (7(1,..., ..., 1) = Ui for all * e {1,..., d}, u t G [0,1], 

(3) For all (ai,..., a d ), (hi,..., b d ) e [0, l] d with a, < h, Ei 1=1 ■ ■ • Z!i d = 1 ( - 1 ) il+ "‘ +t ‘ iC '( u i,* 1 ) • ■ • > u d,i d ) > 0, 
where Uj.± = cij and Ujp = bj for all j £ { 1 ,..., c?}. 

The significance of copulas in studying the joint distribution functions is summarized by the Sklar’s theorem, which shows 
that all joint distribution functions contain copulas and copulas may be used in conjunction with marginal distribution functions 
to construct joint distribution functions. 

Theorem 1. / I29i/ Let F be a joint distribution function with marginals F\,... ,F d . Then there exists a copula C : [0, l] d —>■ 
[ 0 , 1 ] such that, for all x \,..., x d in K. = [—oo, oo] 

F(x u ...,x d ) = C(F 1 (x 1 ),...,F d (x d )). ( 12 ) 


If the marginals are continuous, then C is unique; otherwise C is uniquely determined on RanFi x x ... x RanF d , 

where RanFi = Fi(M) denotes the range of Ft. Conversely, if C is a copula and F[,, F d are univariate distribution 
functions, then the function F is a joint distribution function with marginals F \,..., F d . 


2) Non-Granger Causality: Non-Granger causality is a concept initially introduced in econometrics and refers to a multi¬ 
variate dynamic system in which each variable is determined by its own lagged values and no further information is provided 
by the lagged values of the other variables. This concept has a direct copula expression. 


Assumption 1 . 071 ? The copula function 


C 




■ F 


Xt > • ■ 


. F 


x„ 




■ F 


x„ 


( 13 ) 


admits the hierarchical representation 


( 7*4 (a tj (F fl F t2 F tj d tj (F tr 

L \ X! > r Xi > • • • > XiJ ) • • • > u .Y m \ r X„ 


■ F 


Xrr 


■ F 


X n 


(14) 


where G\. ... ,u{), i = 1 , 2 , ..., m are copula functions. 


Denote the running minimum and maximum of X, up to time tj respectively as 

mifty) = <t k <tj}, (15) 

Mi(tj) = max{Xi(tk);ti <tk <tj} (16) 

and define 

FX(B t ) = P(mi(tj) > Bi), (17) 

F^(Bi) = P(Miftj) < Bi). (18) 

Proposition 1. 071 ? Assume 

PiX^tn) <B 1 , . .., X k (t n ) < B k ) = C tn (F* (Si), ..., F\f k (B k )) , (19) 

and that the copula function allows the hierarchical representation under Assumption [ 7 ] Then the copula function representing 
the dependence structure among the running maximum (minimum) at time t n is the same copula function (survival copula 
function) representing dependence among the levels at the same time, namely, 


P(M 1 (t n )<B u ...,M k (t n )<B k ) = <> (^(S!),...,^^)), 
P{mv(t n ) > B 1 ,...,m k {t n ) > B k ) = C U (f*^ (Bf),..., F^Bkjj . 


( 20 ) 

(21) 
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3) Change of Measure: Consider a stochastic processes {Z t j with a Polish state space E and sample paths in the Skorokhod 
space D = D ([0, oo),E) equipped with the natural filtration {f?t}t> o and the Borel er-field For two processes represented 
by probability measures P, P on ( D ,i^), it is interesting to look for a likelihood ratio process { L t }, such that 

P = E [L t ,A],Ae& t , (22) 

i.e., the restriction of P to (D, j£)) is absolutely continuous w.r.t. the restriction of P to (D, fPf) ll32l l33l . 


Proposition 2. 021/ Let 0 be the natural filtration on D, 3^ the Borel o-field and P a given probability measure on 

(i) If {L t } t >o is a nonnegative martingale w.r.t. ({^t},P) such that E L t = 1, then there exists a unique probability 
measure P on 3P such that ( 1221 ) holds. 

(ii) Conversely, if for some probability measure P and some 3F t -adapted process {L t } t > o ( 1221 ) holds, then {L t } is a 
nonnegative martingale w.r.t. ({J^ t },P) such that E L t = 1. 

Theorem 2. 021/ Let {L t }, P be as in Proposition [2} i). If r is a stopping time and G G 3P T , G C {r < oo}, then 


P(G) = E 


b G 


(23) 


More generally, if the waiting time process W > 0 is fP T -measurable, then E[FF;r < oo] = E [W/L t ;t < oo]. 

Corollary 1. 02l l Let {Lt}, P be as in Proposition}^} i), and let t be a stopping time with P(r < oo) = 1. Then a necessary 
and sufficient condition that E L T = 1 is that P(r < oo) = 1. 

Definition 2. 021/ Assume that {Ztj is Markov w.r.t. the natural filtration 3Ft on D and define {Lt} to be a multiplicative 
functional if {Lt} is adapted to 3Ft and 


Lt+s — Lt • ( L s o Of) 


(24) 


P x -a.s. for all x, s, t, where 6 t is the shift operator. The precise meaning of this is the following: being ,^-measurable, L t 
has the form L t = Tt{{Z u }o<u<t) for some mapping p t : D[ 0, f] [0, oo), and then L s o0 t = ip s ({Z t+u } 0 < u <s)- 

Theorem 3. 021/ Let {Z t } be Markov w.r.t. the natural filtration {^t} on D, let {L t } be a nonnegative martingale with 
E x L t = 1 for all x, t and let P x be the probability measure given by P X (A) = E x [L t \A\. Then the family {P x } x eE defines 
a time-homogeneous Markov process if and only if {L t } is a multiplicative functional. A multiplicative functional {L t } with 
E x L t = 1 for all x, t is a martingale. 

A Markov additive process is defined as a bivariate Markov process {X t } = {( Jt,St )} where {,/ t } is a Markov process 
with state space E and the increments of {S)} are govenrened by { J t } in the sense that 


E [f(S t+s - S t )g{Jt + sWt] = E Jtfi [f(S s )g(J s )\. (25) 

In discrete time, a Markov additive process is specified by the measure-valued matrix (kernel) F(dx) whose ij th element is 
the defective probability distribution 

Fij(dx) = Pi, 0 (Ji = j, Yi e dx), (26) 


where Y n = S n — S n - 1. An alternative description is in terms of the transition matrix P = (pij)ij & E (here pij = Pj(Ji = j )) 
and the probability measures 

Hij (dx) = P(Yi G dx\ J 0 =i,J\= j) = Fi ^ dx \ . (27) 

Pij 

We denote the E x E matrix F[0] with ij\ht element [0] =: J e ex f^A (dx). By Perron-Frobenius theory, the matrix F[0] 
has a positive real eigenvalue with maximal absolute value e K ^ and the corresponding right eigenvector h (0 ' = (h[ S ^) ieE , 
i.e., F[0]h (O ' = e K(B )\f n! . The exponential change of measure corresponding to 6 is then given by 

P = e-^A-^F^Ahts), (28) 

~ e Sx 

H ij{dx) = -£——Hij(dx), (29) 

Hij [u\ 


where A h (s> is the diagonal matrix with the hfb dx on the diagonal, in particular, Pi 3 = e K ^ d )pijh^ fh\, and H i3 [0] is the 
normalizing constant. The likelihood ratio is 


_ b e \j n ) — 6S n -\-nn(6) 

~ h«»(Jo) 


(30) 


which is a mean-one martingale ll32l . 
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III. Copula: A Unifying Technique for Analysis 
A. Bounds of Dependence Structures 

For a random vector X = (Xi,..., X„), where the marginal distribution functions F, ~ X, are known but the dependence 
between the components is unspecified. Denote the sharp lower and upper Frechet bounds over all dependence structures as 
M n (t) := sup {P (]T”=i Xi < t) ; Xi ~ F u 1 < i < n}, m n (f) := inf {P (^,"=1 x i < t )! x i ~ F h 1 < * < n}, and denote 
M+(f) := 1 — m n (t), mf(t) := 1 — M n (t). The sharp Frechet bounds for the case n = 2 were derived in (34) . (35l . and has 
been extended to n > 3, namely the standard bounds Il36l - l38l . 

1) Standard Bounds: 


Theorem 4. 06]? - 0a]? Let Xi ~ Fi, 1 < i < d. Then, for any s G M, we have that 


max 


sup 

i£U(s) 


2=1 


Fi(ui) -(d-l ),0 <P Vli<s < min 


2=1 


inf 

i£U{s) 




(31) 


where U(s) = |u = (m, ..., uf) € : J2i=i u i = s }- 

Remark 2. For the fast computation of standard bounds, a numerical method is described in m while an analytical method 
is described in l\40l . 


2) Dual Bounds: The standard bounds are not sharp for n > 3 and an improvement can be obtained based on duality 
theorems. Specifically, Mf (t) and mf (t) have the following dual counter parts l37l . 


M+(s) = inf <j ^ / gidFt\gi bounded, 1 < i < n with E gi(xi) > l[ s , +oo) E X,; 

V 2=1 ^ *= 1 V 2 — 1 ) 

( n p n /n 

m+ (s) = sup ^ / fidFi ; fi bounded, 1 < i < n with fi(xi) < l[ s , +00 ) ^ X, ; 


(32) 


(33) 


. 2=1 ' 


2=1 


y 2= 1 


While the dual representations are difficult to evaluate in general, they allow to establish good bounds obtained by choosing 
admissible piecewise linear dual functions in the dual problem. 


Theorem 5. £7$, /Ml Let Xi ~ Fj and Fi = 1 — Fi be the survival function of Fi. Then, for any s £ R, we have 

' n r s ~^j^ u i T\l 


M+(s) < D(s) = inf min' 

u€.U{s) 


m n( s ) X d(s) = sup max ■ 

u£U_(s) 


E n rs 

*= 1 Ju. 


Fi(t)dt 


E n 

i= 


, i /, 


i= 1 


n rs-Y. 


E"=t /, 


“3 


Fi(t)dt 


E Tl 

i= 1 U i 


— n + 1 , 0 


(34) 


(35) 


where U(s) ={«£ K n ; ^”_i u i < s l an d U.{s) = {u € R.™; E"=i > s l- 

Let u G U(s), fPfi-i Ui = s, the piecewise linear dual admissible choices become piecewise constant and thus yield a standard 
bound. As a consequence, the dual bounds improve the corresponding standard bounds J37), (38). However, the calculation 
of the dual bounds requires to solve an n-dimensional optimization problem which typically will be possible only for small 
values of n. For the homogeneous case, Fi = F, 1 < i < n, a simplified expression can be obtained with a one-dimensional 
problem that can be solved in any dimension. 

Theorem 6. 071? Let Fi = ... = F n =: F be distribution functions on R+. Then for any s > 0 it holds that 


n f- {n - 1)u F(t)dt 


Mn( s ) ^ D(s) = inf min' 

u<s/n 


rrif ( s ) > d(s) = sup max < 

it>s/n 




s — nu 

s—(n— l)u 


F(t)dt 


,1 


— n + 1,0 , 


s — nu 

Theorem 7. 44/1 ? Let F\ = ... = F n =: F be distribution functions on M + . Then for any s > 0 it holds that 

ix 


(36) 

(37) 


m n (t ) > 1 — n inf 

r£[0,s/n) 


f r s - (n ~ 1)r F( X )da 


s — nr 


(38) 
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The infimum in ( |38| > can be easily calculated numerically by finding the zero derivative points of its argument in the specified 
interval. Note that 

)r F(x)dx' 


lim < 1 — n- 

r—fs/n I 


= nF(s/n) — n + 1, 


(39) 


which means that (f38l is greater or equal than the standard lower bound. Moreover, the dual bounds in (l38l > have been proved 
to be sharp under some general distributional assumptions l42l : 

(Al) Attainment condition: There exists some a < s/n such that 


D(s ) = inf 

t<s/n 


n Jt 


-7=7 


F(x)dx n F(x)dx 


s — nt b — a 

where b = s — (n — 1 )a and a* = F _1 (l — D(s)) < a. 

(A2) Mixability condition: The conditional distribution of ( Xi\Xi > a*) is n-mixable on (a, b). 
(A3) Ordering condition: For all y > b it holds that 


(40) 


(n-l)(F(y)-F(b))<F(a)-F 


s-y 

77—1 


(41) 


Theorem 8 . P2l/ Under the attainment condition (Al), the mixing condition (A2) and the ordering condition (A3), the dual 
bound is sharp, that is 

n It ^ 1 '* t F(x)dx nJ^F(x)dx 


M+{s) = D(s)= inf —-------. 

t<s/n s — nt b — a 

The first order conditions for the optimization in the attainment assumption (Al) at t = a imply that 

rb~Fi, 


n f a F(x)dx 


b-, 


= F(a) + ( 77 - 1)F(6), 


(42) 


(43) 


where b = s — (n — l)a and a* = F X (1 — D(s)) < a, and it provides a clue to calculate the basic point a and, hence, the 
dual bound D(s). Having calculated a one can easily check the second order condition 

f(a)-(n-l) 2 f(b)>0, (44) 


which is necessary to guarantee that a is a point of minimum for (Al). At this point, the sharpness of the dual bound D(s) 
can be obtained from a different set of assumptions l42l : 

1) Continuous distribution functions F having a positive and decreasing density / on (a*,oo) satisfy assumptions (A2) and 
(A3). 

2) Continuous distribution functions F having a concave density / on the interval (a, b) satisfy the mixing assumption (A2). 
In order to obtain sharpness of the dual bound D(s) for these distributions, conditions (Al) and (A3) have to be checked 
numerically. 


B. Example Dependence Structures 
1) Comonotonicity: 

Definition 3. H43\l The set A (— I9L 11 is said to be comonotonic if for any x_ f /y or y f x_ holds, where x_ y denotes the 
componentwise order, i.e., Xi < yi for all i = 1,2,... ,n. 

Definition 4. fi43\l A random vector X_ = (X\,..., X n ) is said to be comonotonic it has a comonotonic support. 

From the definition, we can conclude that comonotonicity is a very strong positive dependency structure. Indeed, if x and 
y are elements of the (comonotonic) support of X_, i.e., x and y are possible outcomes of X, then they must be ordered 
componentwise. This explains why the term comonotonic (common monotonic) is used l43l . 

Theorem 9. f43d A random vector X_ = (X\, X 2 , ■ ■ ■, X n ) is comonotonic if and only if one of the following equivalent 
conditions holds: 

(1) X_ has a comonotonic support. 

(2) For all x = (xi,X 2 , ■ ■ ■, x n ), we have 

Fjcfe) = min{F Yl (ar), F X2 (x 2 ),..., F Xn (x n )}. (45) 

(3) For U ~ Uniforin(0, 1), we have 

X ± (F^(U), F ^([/),..., F x *(U)). (46) 

(4) There exist a random variable Z and non-decreasing functions fi(i = 1,2,..., n), such that 

X±(f 1 (Z),f 2 (Z),...,f n (Z)). 


(47) 
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2) Independence: 


Theorem 10. jf29|7, fJUf Random variables with continuous distributions are independent if and only if their dependence 

structure is given by 

d 

C(ui,...,u d ) =Y\_Ui- (48) 

i= 1 

3) Markovian: The Markov property is a pure dependence property that can be formulated exclusively in terms of copulas, 
as a consequence, starting with a Markov process, a multitude of other Markov processes can be constructed by just modifying 
the marginal distributions mi, mu. 


Definition 5. [WH, fcgH, fi4$l Assume two bivariate copulas A(u,t) and B(t,v ): the product operator * is defined as 

nl dA{u, t) dB(t, v) 


A * B(u, v) — ( 

Jo 


-dt. 


dt dt 

Assume an m-dimensional copula A and an n-dimensional copula B: the start operator * is defined as 

dA(u \,.. 

• i 'U'm—l ? dB(t , 

5 ■ ■ • 5 tlm-yn— 1 ) 


A * B(lli , U2 , • ■ ■ , Um-j-n—i ) — 


dt 


dt 


dt. 


(49) 


(50) 


Theorem 11. fWj, §44%, /g5]? A real valued stochastic process X t is a Markov process of first order if and only if for all 
positive integers n and for all ti, ... ,t n satisfying tk < tfe+i, k = 1 ,..., n — 1 , 


C t - 


,,t n — Cti,t 2 Ct 2 ,t 3 * . . . * Ct n _ 1 ^t ri 


(51) 


where is the copula of X tl ,... ,X tn and Ct k ,t k+1 is the copula of X tk and X tk+1 . 

A strictly stationary Markov process is characterized by assuming C\ i _ 1 \ x = C, Vi and F\ = F, Vi. This kind of process 

dC , 
du 


can be constructed for any given C and F, with D±C(u,v) = tFf(u,v), by setting 


and 


Fxi(t) = f D\C{w, F(t + F 1 (w)))dw = G(t), 

Jo 

Cx^Xi («,«)= [ DiC(w, F(G _1 (t) + F~ 1 (w)))dw. 
Jo 


Definition 6. fl7\l Assume that, for x £ [0, l] fc the distribution C(x,.) is absolutely continuous with respect to the 
generated by some copula A on [0,1] ( . We denote by Ca(x, y) (a version of) theRadon-Nikodym derivative, 

C( x , dy) = C iA (x,y)A(dy). 


(52) 

(53) 

isure 

(54) 


The subscript “, A” indicates that we take the derivative with respect to the second set of arguments (y). Accordingly, we 
define the derivative Cb,(x, y) of C with respect to a k-dimensional copula B by 


C(dx., y) = C B ,{x,y)B(dx), 


(55) 


provided that for given y the measure generated by C(., y) is absolutely continuous with respect to the measure generated by 
B. Cxjxj) and Cb,(x, y) are called derivative of the copula C(x, .) resp. C(., y) with respect to the copula A resp. B. 

Definition 7. 4771/ Let A be a (k + m)-dimensional copula, B be a (m + l)-dimensional copula and C be a m-dimensional 

C(.) 

copula such that the derivatives A : c and Be, are well-defined. The operator ★ is defined by 

(A C * ) B)(x,y) = f Ac(x,r) • B c ,(r,y)C(dr), (56) 

provided that the integral exists for all x, y, z. 

Theorem 12. The n-dimensional process X is a Markov process if and only if for all t\ < t 2 <■■■< t p the copula 

c t u ...,t p (x tl ,..., X tp ) satisfies 


CT 


, t t C*2(.) t . C* 3(.) 

;t p=C t 1 ,t 2 2,t 3 f-’ 






- 1 J tp 


(57) 


Given a Markov family C st , s.t G AT, s < t, we can define finite dimensional copulas C tl ’"" tp and combine these with 
an arbitrarily specified flow F t (x) = (F x i(xi),..., Fxy (x n )) of marginal one-dimensional distributions F X i to obtain finite 
dimensional distributions 


P(X tl <x 1 ,...,X tp <x p ) (F tl ( Xl ),...,F tp (x p )). 


(58) 
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By applying Kolomogorov’s construction theorem for stochastic processes there exists a Markov process X with the given 
copulas and marginal distributions. 

Remark 3. The copula representation of Markov processes of order k has been extended in mi . mi and the construction 
of Markov processes through increment aggregation is elaborated in l\45d . Moreover, the results for the one-dimensional case 
in nBi are extended to the general multivariate setting in 4771/ . In addition, the relationship between Markov property and 
martingale is investigated in m-m. 

Theorem 13. f44\1 f45\l Let X = {Xj}j> 0 be a Markov process and set Yi = Xi — Xj_i. X is a martingale if and only if: 

(1) F Yi has finite mean for every i; 

(2) for i > 1, Fy. 1 (v)dD 1 Cx i _ 1 ,Y i (V',v) = 0, Vu <E [0,1]. 

Proposition 3. ff4h . Jf45j/ Any process for which the distribution of increments is symmetric and the copula between the 
increments and the levels is symmetric (around the first coordinate) is a martingale. 

By definition, (X, Z) is a martingale with respect to JF X > Z iff Vi > 0 


E[Xi+!-XJX^Zi] = 0, (59) 

— Zi\Xi, Zi\ = 0. (60) 


Let A Xi = Xi + i~Xi, AZi = Zi + i—Zi and A li+ i ( u, v, w, A) be the copula function of the random vector (Xi, Zi , A Xi, AZf) 
with Fxi , Fz t , F/sXi, F AZt the corresponding marginal cdf. Set a,iy+i(u, v, w, 1) the density of the copula Ai^ + i(u,v,w, 1) 
and v, 1, w) the density of the copula A+i+i(it,i>, l,w). 


Theorem 14. f45\l The bivariate Markov process (X,Z) is a martingale with respect to the filtration & x,z iff: 

(1) F/sXi and Fxz, have finite mean for every i; 

(2) for every i, \/u,v £ [ 0 , 1 ], 


F AX t ( w ) a i,i +1 ( u , V, W, 1 )dw 

= 0, 

(61) 

F AZi ( w ) a i,i +1 ( U , V, 1, w)dw 

= 0. 

(62) 


IV. Analysis of Cumulative Capacity 


A. General Results 


1) Exact Expression: The CDF of the cumulative capacity can be derived from the joint distribution of the channel gain 
magnitude, i.e., 


Fs(s,t)(x) = l 


fn(h s +i, h s+ 2 , ■ ■ ■, h t )dH. (63) 

/S(s,t)=X)* =s+1 log 2 (l+7|ft,i| 2 )<a; 

For example, the single-integral form PDF and CDF of the multivariate generalized Rician distribution are expressed as 


N 


/h((ii, ( 12 , ■ ■ ■ , hx ) = 


— exp(-(f + S 2 ))I m -i(2SVt) „ t 


'4=0 


1 


X ^2 /l ™ eX P I “ 


hl + xMt 


2 


dm— 1 I dk 


=1 iKo-l-tm 
y/v'iKt' 




dt, 


Fa(hi, /i 2 , • • •, /ijv) = 


m — 1 

t 2 


N 


't=0 


— exp(-(t + S 2 ))I m -i(2SVt) U 


fc=1 L 


1 ~Qr 


v / ^V / °'fc A fc hk 


dt, 


(64) 

(65) 


where = " 2 ‘ 1 A 


According to 


and S 2 = ( m h + m li)- 

the transform of random vector that is 

F s(s,t)( x ) = / 

Js 


is elaborated in Theorem [15] the CDF can also be expressed as 

/c(t's+i? C3+2? ■ ■ ■ 1 Ct)dC, 


( 66 ) 


' S(s,t)=J2i= a +l c i< x 
where / c (y) = /h(C'” 1 (y))| J 0 (y)| and C(y) = log 2 (l + 7 |H(y)| 2 ). 

Theorem 15. / |49l) Let g : M" -A R™ be one-to-one and assume that h = g 1 is continuous. Assume that on an open set 
V C R” h is continuously differentiable with Jacobian J(y). Define J 0 '■ R n —> R" by 

J(y), 

0, y£V c 


Jo{y ) = 


(67) 
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where V c is the set complement (in R") of V. Suppose random vector X has pdf fx(x) (with respect to Lebesgue measure) 
with nonzero mass in h(V c ), i.e., P{X £ fr(V c )} = f vc fx{x)dx = 0. Then the pdf of Y = g(X) is given by 

/,(») =/x( s -'fa))l*(ril = (^ (r,(! ' ))|J<! ' )1 ' m 

[0, y£V c . 

In addition, the exact formula can also be expressed with copula. For a d-dimensional random vector X, denote pc the 
measure on [0, l] d implied by the copula C, that is pc{B) := P[U £ B] for any Borel measurable B C [0, l] d and U ~ C. 
Define Z = X]f=i <^iXi, let Az be the convex set Az = {x £ R d : u i x i < z }, and A* z be the linear boundary connected 

to A z , i.e., A* z = {x £ R d : Yh=i ^iXi = z}. 

Proposition 4. K50\l 

F z (z) = pc(Bz) = [ cdX , (69) 

J Bz 

where c is the density of C, Bz = <p(Az) C [0, l] d , and <p : R d —>• [0, l] d , (x\, ..., x d ) H > (Fi{xi), ■.., Fd(xd))- 

Proposition 5. K5W Assume that f , > 0 on R for i, = 1..... d, then 

(1) the set B* z = ^>(A* Z ) is given by B* z = {u £ [0, l] d : (tti,.. .,u d - i) £ (0, l) d_1 ,u d = t z (u i, ... ,u d -\)}, where 
r z : (0, l) d_1 -f (0, l),u£ F d (^ d - E-=1 ^ _1 K)); 

(2) the set B z is given as B z = {u £ [0, l] d : (m,... ,u d - 1 ) £ (0, l) d_1 ,0 <u d < t z (ui, ... ,u d - 1 )}; 

(3) for z i < z% B Zl C B Z2 ; 

(4) r z (u) is a strictly decreasing function in each component of u; 

(5) B z is path-connected and B z := B z D [e, 1 — e] d is start-shaped for e > 0 with center (e,..., e). 

Remark 4. Algorithms for numerical calculation of Fz(z) are elaborated in ( I5QI/ - I52]/ . 


2) Bounds: In the following, instead of trying to find accurate representation of Fs( s ,t), we turn to investigate bounds on 
it. With Theorem |4] it is easily verified that the CDF of the channel cumulative capacity satisfies the following inequalities: 


F l s(S:t) (r) < F s{Sit) (r ) < F^ s t) (r), 


(70) 


where 


= t inf 
Z) r *= r 

i = s + l 


F S(s,t) ( r ) = Slip 

t 

Z) r i= r 

i = s + l 


Z F c(i)(n ) , 

_i=s+l J i 

" t 

Z F C(i)( r i ) -(t-s- 1) 

_i=s- 1-1 


(71) 


(72) 


The improved results with dual bounds can also been obtained following Theorem [5] to Theorem [8] 


B. Special Cases 

1) Comonotonicity: The copula of a distribution function contains all the dependence information, specifically for a 
comonotonic dependence structure, i.e., comonotonic random variables are increasing functions of a common random variable 
||39l , the joint distribution function of the instantaneous capacity is expressed as ED, ED 

F(C(s + 1), ..., C(t)) = min {F(C( S + 1)),..., F(C(t))} , (73) 

and the CDF of the cumulative capacity is expressed as ED, ED 

m t -s{ x) = [ dmm{F(C(s + l)),...,F(C(t))}. (74) 

JY.\=s +i C(i)<x 

In the special case that all marginal distribution functions are identical -Fc(i) ~ Fq, comonotonicity of C(i) is equivalent 
to saying that C(s + 1) = C(s + 2),..., = C(t) holds almost surely ED, in other words, the capacity depends only on the 
initial state and keeps constant afterward, i.e., 

F S (,,t){x) = F C (^)- ( 75 ) 
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2) Independence: If C(i) and C(j ), i ^ j, are independent, fs( s ,t) = fc(s+ 1) * ■ ■ • * fc(t)i where * denotes the convolution 
operation. Hence, 

Fs(s,t){x ) = f fs(s,t)(y)dy- (76) 

J —oo 

In addition, when the length of the period, t — s, is large, according to the central limit theorem, I's(s.t) ( x ) approaches a 
normal distribution with mean £J[5(s,f)] and variance cr 2 [>S'(s, f)](under certain general conditions) 1531 . i.e.. 


Fs(s,t)(x) ~ G ^ 


x — E[S(s, f)] 
cr 2 [5'(s,f)] 


(77) 


where E[S(s,t)] = £ E [ c< d)i cr 2 [5'(s, t)] = cr 2 [C(i)], and G{x) = e y2 ^dy. 

i=s+l 2=S+1 

In the special case that all marginal distribution functions are identical Feu) ~ Fc, the cumulant generating function and 
the likelihood ratio are expressed as 


k( 0) = logEe ec(i) =log J e ex F(dx), 

T _ 6St — tn(6) 

■ Lj t 5 

where L t is a mean-one martingale. According to the Markov inequality 


P{L t > p} < — E[L t ] — —, 


i.e., for 9 < 0 


P{S t >x}< e ex - tK ^. 

3) Markovian: The likelihood ratio martingale can be used to provide exponential upper bounds f54fl . Define 


= 


min n {h^ e \j n )) _ 


0Sn~\-nK(0 ) 


(78) 

(79) 

(80) 

(81) 

(82) 


hW(J 0 ) 

Then L n < L n , so E[L n ] < 1. For any p > 0, the Markov inequality and Doob’s submartingale inequality can be applied to 
L n and L n to get 


P{L n > p} < i]E[Lj < 

/i /i 

P i sup L n > p\ < -E[L n ] < 
U>i J p p 

For any a and 9 < 0, we can choose p = e nK ^~ ea mm n (h^(J n ))/h^ e \Jo) 


P{S n > a}< 


h^(Jo) 


min n {h^\j n )) 


o —nK,(0)-\-0oi 


(83) 

(84) 

(85) 


V. Analysis of Maximum and Minimum Cumulative Capacity 

A. General Results 

With assumption of non-Granger causality on the dependence structure, a lower bound for the CDF of the maximum 
cumulative capacity and an upper bound for the CDF of the minimum cumulative capacity are obtained 


sup S(i) < x 

0 <i<t 


P (5(1) < x, 5(2) S(t) < x) 


> 


oc 

max C(i) < — 
1 <i<t t 


c (f m 1 {x),Fm 2 ,..., F Mt ) 

c(r W ,r(|.|),...,r(”,...,£)). 


( 86 ) 

(87) 

( 88 ) 

(89) 





P ( inf 

\0<i<t 


S{i) < 





12 

P(S(1) > x,S( 2) > x,. 

. .,S(t) > x) 


(90) 

P ( minC(l) > x, min 

V l<i<2 

C(i)> 

min C(i) > — ) 
l<i<t t J 

(91) 

c(F mi (x),F m2 (|), 

(j)) 


(92) 

c(T,),r (If),... 


?))■ 

(93) 


where F{x 1 ,x 2 ,...,x t ) = C(F c{1) (xi), F c(2) (x 2 ), ■ ■ ■, F c{t) (x t )). 


B. Special Cases 

1) Independence: In the special case that all marginal distribution functions are identical Fqu) ~ Fq, the cumulant 
generating function and the likelihood ratio are expressed as lp5 " =rl 


re(0) = logEe 00 ^ = log J e ex F(dx), (94) 

L t = e eSt ~ tKie \ (95) 

where L t is a mean-one martingale. Let the Lundberg equation k(0) = 0 and assume the existence of a solution 9 > 0, then 


P <! sup S t > x }> < e~ 9x 
,t> o 


(96) 


for all x > 0. 


Jl U.X1 Oj ^ . 

2) Markovian: Let t(u) = inf{f > 0 : S t > u }, I(u) = J T ( U ), £(u) = S T ( U \ — u, M = sup t>0 S t . Let the Lundberg 
equation k(6) = 0 and assume the existence of a solution 9 > 0. Then ir ^ 5 ’ 1 


P i{M > u) = P i(r(u) < oo) = E it e 


h {6) 

n J 0 „-es T(u) 


.(B) 

2 Je(.u) 


; t(u) < oo 


h) 


(B) 


h (B) 

% I(u) 


= e~ 9u ^ e 
P (M > u) = TtjFj. 

i 

According to Lundberg’s inequility l55l 

P i(M > u) < 




l(«) 


„— 6 u 


min j eE hj 
P [M > u) < 7TiPj. 


(«) 


(97) 

(98) 

(99) 

( 100 ) 

( 101 ) 


The constant in front of e 9u in Lundbergs inequality can be improved and a supplementary lower bound can also be 
obtained with the following theorem. 

Theorem 16. j55l/ Let 


1 . f Bj(x) 

C- = min-Tjr ■ ml ——-— - 

3&e p(B) x>o f™ e e (y~ x )Bj(dy) 


C+ = max 


1 


Bj{x) 


= max—r- • sup rOC .. ~ —r ——, 

h (B) x >o j™ e^v-^B^dyY 


where Bj is the distribution of the increment. Then for all j £ E and all u > 0, 

C-h\ e) e~ eu < P i(M >u)< C+h[ e) e~ eu . 


( 102 ) 

(103) 

(104) 


Remark 5. The above results in this section have been so far essentially for the forward-looking maximum / minimum 
cumulative capacity. For their backward-looking counterparts, similar analysis may apply, and so may for the the general 
definitions of maximum and minimum cumulative capacity. 
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VI. Other Characteristics of Cumulative Capacity 
A. Moment Generating Function 

t 

The moment generating function (MGF) of S(s,t) = (7(f), denoted by Mg( st )(0), is (5), |f71 

i=s+l 


M S{s ,t){0)=E 


0 es(s,t) 


f 


e 0r dF S{Sit) (r), 


(105) 


where 0 is a real variable. We denote = Afs(«,t)( — 0) = -E [e -6,s ( s ’0]. 

With F S ( S: t) (x) derived in the previous section, Ms( s ,t)(&) I s readily obtained. Specifically, for the independent case, we 
have 



M l s ^ t) {6) = nU +1 M c{i) (0) 

(106) 


„ r e -«-dG( r - E[SM] '\ 

J-oo l ^[S(s,t)} )■ 

(107) 

where for the independent and identically distributed (i.i.d.) case, (|l06|l becomes 



M il SM {6) = (M CW (0))‘- S . 

(108) 

For the more general case. 

Mf (Sit) (0)<M^ jt) (0)<M^ M) (0), 

(109) 

with 

hi r°° 

M S{ s,t)(0) = J J~ er dF l SM {r\ 

(110) 


= r e--d^ (s , t) (r). 

J — OO 

(111) 


1) Rayleigh channel as an example: For the i.i.d. case. 


M’s(.,t)(f)= (j[°°e- flr d{ j™ e~ er 2 r e-^-^dr ^' 


( 112 ) 


For the general (possibly dependent) case, the upper bound and the lower bound of the MGF M S ^ s t )(0) are expressed as 

-f 2 *— l) /7 


—du 


MslM 9 ) = / e_Srd 

Jo 

= (t — s) ( 1 — exp 


(f — s) ^1 — e 

X _ 2 n “/0 _s ) 


— Ofh 


7 


— 0(f — s) exp 


1 _ 2 r /( t-s ) 


7 


— Or I dr, 


-di 


M s\s,t)( e ) = J Q e erd ^ 1 -(f-s)e 

'1 _ 2 n i/( t “ s ) 


= (f — s) exp 


7 


-( 2 ‘- s - 1 )/ 

— — 0(f — s) J exp ^ 


1 _ 2 r /( t_s ) 
7 


— Or I dr. 


(113) 

(114) 

(115) 

(116) 


where il u = rlog 2 (1 — 7 log(l — 1/r)) for r > 1, fl; = r log 2 (l — 7 log(l/r)) for r > 1, and r = t — s. 

2) Application to effective capacity: Suppose the service process has stationary increments, i.e., S(s,s + t) = st S(t) for all 
s,t < 0, the effective capacity of S is defined as J6|, l56l 


1 


F c) (0) = -lim sup*^ — sup log E 
"t s >0 


0 —B(S(s,s+t)) 


where 0 is a real variable. 

The effective capacity calculated through the instantaneous capacity follows directly from (1 1 07b and (1 1 09b . i,e., 


= — limsupt^-j-suplog Mj (SiS+t) (0), 


1 


Of s>0 


(117) 


(118) 


r { d\Q) = -limsup^oo^-suplog M^ Sig+t) (0), 


Of S>0 


(119) 
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for the independent and dependent cases respectively. 

For the Rayleigh fading channel, log Mj( T ) ( 9 ) is constant with r, ^ log M s \‘ T ^ ( 9 ) increases with r, while log M S( - r ) ( 9 ) 

decreases with r. The effective capacity for the i.i.d. case and dependent case are respectively, 

r&°(0) = -- Q log J e -er 2 r e -( 2r - 1 )/ 7dr ^ ; (120) 

r d ?( 9 ) = - ^lim^ y log (V(l - exp((l - 2 n “ /T )/7 - 9 £l u )) - 9t f exp ((1 - 2 r/r )/7 - 9 r) dr) , (121) 

J 0 

r { d u( 9 ) = -ilog(rexp((l-2 n! / T )/7-0F! i ) exp ((1 - 2 r / T ) 7 - Or) dr). (122) 


B. Mellin Transform 

If the channel capacity is expressed with the natural logarithm, i.e., 

C(t) = log(5(7*)) nats/s/Hz, 


(123) 


in order to circumvent the calculation of logarithm, the cumulative capacity may be transformed to a new domain, e.g., the 
MGF domain, with the exponential. 


S(s,t)=e s ^= 5 ( 7 *)' 


(124) 


i=s+l 


Then the MGF property analysis applies. 

Analog to the MGF for the cumulative capacity, Mellin transform (MT) has also been proposed fT 5 l , i.e., 


Af s(Sit) (d) = 

/ oo poo 

r*- 1 dF s[s , t) (r)= / 

-OO J —OO 

where r* = log 2 e log r, for any complex variable r) for which the right hand side of 

According to ( 1771 ) and (fTOl . the MT of the cumulative capacity for the independent and dependent case can be derived, 

1^ (r*-E[S(s,t)] 


(125) 

(126) 


exits. 




5(s,t) W 


--o/" dG \ <? 2 [S(s,t)} 




< 


dF l SM (r*) < M d S(S ' t) (#) 


S(s,t ) \ 


/ oo 

-OO 


where ( 11281 ) and ( 11291 ) hold when 1? — 1 < 0 . 

1) Rayleigh channel as an example: For the i.i.d. case tm 




For the dependent case. 


Mill'd) = I r^d- 


log 2 log 2 e 


(f — s) ^1 — e 

„ P n u/ l°g2 e 


- 2*— -1 / 7 


tf -2 


x 2‘-= x e t dr , 


(127) 

(128) 

(129) 

(130) 

(131) 

(132) 


r‘ d ~ 1 d • 


1 — (t — s)e 


2‘—-1 J / 7 


>J rxe^~ 2 )r d' 


1 — (f — s)e 


- 2*-*-i )h 


log 2 log 2 e 
7 


x / e' 

e p. l /\og 2 e 


o r)\ r* 1 —2 * — s 

? -A x 2*- s x e -v dr, 


(133) 

(134) 

(135) 


where d — 1 < 0, f2„ = rlog 2 (1 — 7log(l — 1/r)) for r > 1, 12; = Tlog 2 (l — 7log(l/r)) for r > 1, and r = t — s. 
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C. Stochastic Strict Service Curve 

A system is said to be a stochastic strict server providing stochastic strict service curve (SSSC) P with bounding function 
g, if during any period (s, i] the amount of service S(s,t) provided by the system satisfies [0, 

P {S(s, t) < f3(t — s) — x} < g(x), (136) 


for all x > 0. 

1) Independent Case: If C(i) and C{j), i / j, are independent, applying (f77l) . the expression of SSSC is obtained 

' iP{t - s) - x) - E[S(s,t)]' 


P{Sis,t)<p(t-s)-x}nG 


= 9lA x )- 


a 2 [S{s,t)\ 

If S(s,t ) also has stationary increments, i.e., C(i) = s t C(j), the SSSC can be expressed as |8) 

P l {t- s) = ^log M l S{Stt )(9) = t -^- log M c (6), 


with bounding function g(x) = e . 

2) Dependent Case: If S(s,t ) has stationary increments, the SSSC can be expressed as (SJ 


P d (t-s) = L log M d s{S ' t) (8), 


(137) 


(138) 


(139) 


with bounding function g(x) = e~ 6x . 

More generally, if C{i ) and C(j), i / j, are dependent, replacing r = p(t - s) x in (1711 and (l72l) . the SSSC of the 
capacity is obtained 

= F s( s ,t)(P(t -s)-x) <P{S(s,t ) < Pit - s) - x} < Fs( s , t )(P(t - s) - x) = gfA x )- (140) 

3) Rayleigh channel as an example: For the i.i.d. case, applying (11381) . 


P U (t -*) = *—£ log Mc(0) = ^ log (jj- J™ e- 0r 2 r e-( 2r ~ 1 ^dr) , 

with bounding function g(x ) = e~ 8x . 

If S(s,t) has stationary increments, the lower bound and upper bound of the SSSC can be expressed as 


P dl (r) = -^log M? (r) (0) 


— log (r(l - exp((l - 2°“ /t )/ 7 - 0fl u )) ~ dr J ex P ((! ~ 2? " ,T )h ~ dr) , 


1 


-dl 


P du ( T ) = ^log M s[t) { 6) 


-9 
1 


L log (r exp ((1 - 2 n ‘/ T )/ 7 - 0fi,) - 9t £ exp ((1 - 2 r / T )/ 1 - Or) dr) , 


where t = t — s, with bounding function g(x) = e~ 6x . 

For the dependent case, the lower bound and upper bound of the SSSC are expressed as 

t 


P {S(s,t) </3(t — s) — x} < inf 


E ri=P(t-s)-x 
= 0 + 1 




,i=s+l 


< 


(t~s) (l-exp(- (Wi -l)/ 7 ] 


■))] i =5^( a; ), 


P {5(s, t ) < P{t — s) — x} > sup 


E n 

>=o+l 

= /3(t — s) — x 


E (l-e-(» r ‘-«/7) 1) 


,2 = S +1 


> 


1- (t- 


it — s) exp ^2 


/3(t — s) — x 


-i h 


i + 




Let g1: u t (x) = e, r = t — s, and Pf{r) = j3(r) — x , we obtain 


Let gf t (x) = e, r = t — s, and /3^(r) = /3(r) — x, we obtain 


Pl (r) = t log 2 ^ 1 7 log (l- ^)) ■ 

r) — x , we obtain 
Pi(r) =r\og 2 fl- 7 log 


1 — £ 


(141) 

(142) 

(143) 

(144) 

(145) 


(146) 

(147) 

(148) 

(149) 

(150) 

(151) 
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VII. Conclusion and Future Work 

The exploding data and increasing integration in future wireless communication put forth the exploration of the ultimate 
capacity that the wireless channel can provide while ensuring more and more stringent demand on quality of service. 
Novel information theoretic measures other than the conventional ergodic capacity or instantaneous capacity are required 
for network performance analysis. This paper has advocated a set of (new) concepts to study or as further properties of 
wireless channel capacity. They include cumulative capacity, maximum cumulative capacity, minimum cumulative capacity, 
and range of cumulative capacity, where for the latter three concepts, their forward-looking and backward-looking variations are 
also discussed. In addition, copula has been introduced as a technique to unify analysis of these introduced channel capacities 
under different dependence structures of the underlying cumulative channel capacity processes. Extensive results, including 
exact solutions for special cases and bounds for more general cases, for cumulative capacity are derived, mainly in terms of 
probability distribution function characteristics, e.g. CDF of the cumulative capacity process. In addition, its other characteristics 
in terms of moment generating function, Mellin transform, and stochastic service curve are also investigated. For maximum 
cumulative capacity and minimum cumulative capacity, preliminary study and results are presented. 

It is an on-going work is to extend the study on maximum and minimum cumulative capacity and use the range of cumulative 
capacity as a measure of tightness of the upper and lower bounds of the cumulative capacity. Both forward-looking and 
backward-looking variations of the general definitions are being studied. This work will be extended to apply the introduced 
concepts to QoS performance analysis of wireless channels, in addition to QoS-constrained capacities e.g. delay constrained 
capacity (57H59). At present, the analysis in this paper is based on a default assumption that the dependence structure of the 
cumulative service process is stationary over time. If the dependence behavior is time-variant, dynamic copula should be used 
to used for dependence modeling lf60l . which will be our future work. 
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